Full counting statistics of Kondo-type tunneling in a quantum dot: the fluctuation 

eff'ect of Slave-Boson field 
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We study the full counting statistics (FCS) of electron tunneling through a multi-terminal quan- 
tum dot in the Kondo regime within the slave-boson mean field theory. By employing the A.O. 
Gogolin and A. Komnik's method of calculating the FCS generating function based on the nonequi- 
librium Green's function [Phys. Rev. B 73, 195301 (2006)], we obtain the counting field A-dependent 
self-consistent equations for the mean values of the slave-boson fields and the explicit expression for 
the derivative of the adiabatic potential of the system with respect to the counting fields. Perform- 
ing perturbative expansion to the first order of A, we find an extra contribution to the shot noise 
due to the bias-induced Bose field fluctuation, and then confirm that the nonequilibrium particle 
number fluctuation plays an important role in the current noise of the Kondo dot: enhancement of 
the current auto-correlation and a positive current cross-correlation. 

PACS numbers; 72.10.Fk, 72.70.-|-m, 73.23.-b, 73.63.Kv 
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I. INTRODUCTION 

Recently, the full counting statistics (FCS)ii^ of charge 
transport in mesoscopic system has become an active 
topic of experimental^Ti^ and theoretical^Ti^i investiga- 
tion. This concept was first proposed by Levitov and 
Lesovik^ to describe the whole probability distribution of 
transmitted charge during a fixed time interval in a non- 
interacting mesoscopic conductor. It is therefore beheved 
that the FCS contain full information about electron cor- 
relations that can not be obtained via measuring the av- 
erage current alone^ii^ Since then FCS has been studied 
using the scattering matrix theory in a variety of systems, 
for example, normal and superconductor hybrid tunnel- 
ing junctions ^SrJ^ chaotic quantum dots (QDs),^^ entan- 
gled electrons^ and spin current fiS, Moreover, for taking 
account of the effect of electron-electron interaction, a 
general scheme has been developed to evaluate the FCS 
of the Coulomb-blockade QD in the framework of a quan- 
tum master equation approach. This scheme has been 
then applied to analyze the FCS in a multilevel QD^ii 
a coupled QDs,— '-^ and a single-molecule magnet. In 
combination with standard Rayleigh-Schrodinger pertur- 
bation theory, this approach has been further improved 
for investigation of the FCS in nano-electro-mechanical 
systems2i and even for consideration of non-Markovian 
dynamics.— Very recently, finite-frequency FCS has been 
explored in an interacting QD.^^ Besides, one of the 
authors has extended the famous MacDonald's formula 
to calculate the FCS for the nonequilibrated-vibration- 
assisted tunneling in a molecular QD.^^ 

Over last decade, since its experimental discovery in 
nanoscale devices the Kondo effect and related physics 
have attracted enduring attention. This is owing to the 
impressive advantage of an artificial atom, easy tunabil- 
ity of system parameters in a considerably wide range and 
full controllability of external circumstances, which facil- 
itate a new realm of research, the nonequilibrium Kondo 
physics in transport measurements ?SS^^^ For instance, a 



bias- voltage-induced splitting of the Kondo peak in the 
local density of states of the QD has been experimen- 
tally observed in a three-terminal transport setup. It 
is therefore desirable to explore the FCS of electron tun- 
neling passing through a QD connected to multi-terminal 
in the Kondo regime. However, those approaches men- 
tioned above are not applicable for this purpose. 

Actually, many attempts have been made in the lit- 
erature to study the zero-frequency shot noise, i.e. the 
second current cumulant of the FCS, in a Kondo-QD<2^ 
Hershfield^- calculated the zero-frequency current noise 
using perturbation theory in the Hartree approximation 
based on the Green function (GF) approach. He found 
the interaction can either enhance or reduce the shot 
noise. Yamaguchi and Kawamura'^*' performed a comple- 
mentary analysis by choosing the tunneling term of the 
Hamiltonian as the perturbation parameter and revealed 
a large suppression as comparison with the Poisson value. 
Ding and Ng^i calculated the frequency-dependent shot 
noise by employing the equation-of-motion method and 
Ng's ansatz for the correlation GF. Their results also 
demonstrated suppression of the shot noise below the 
noninteracting value. However these results are all only 
valid at high temperature, T > Tk (Tk is the Kondo tem- 
perature), and weak Kondo correlations. For analyzing 
the shot noise in the case of strong Kondo correlation, 
Dong and Lei^ employed the finite Coulomb interaction 
U slave-boson mean-field theory (SBMFT) to calculate 
the current-current correlation function at zero temper- 
ature based on the nonequilibrium GF (NGF) technique 
and Wick's theorem. Later, Lopez and his coworkers uti- 
lized the same theoretical framework with the infinite-C/ 
version of SBMFT to study the shot noise in a single 
QD with ferromagnetic leads^^ and multi-terminals?^ 
and a coupled QD system.^ Very recently, this method 
has been used to evaluate the shot noise in the parallel 
coupled QDs to distinguish between the spin and orbital 
Kondo effect 1^ It should be noted that the fluctuations 
of slave Bose fields are all neglected in the calculation 
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of current correlation function. 32, Meir and Golub per- 
formed a exhaustive calculation of the shot noised based 
on the noncrossing approximation (NCA). However, they 
just substituted the resulted NCA propagators into the 
current correlation function derived from the noninter- 
acting electron presumption. 

For applying the SBMFT to explore the FCS, one 
must overcome a technical difficulty: how to evaluate all 
higher-order terms of current correlation functions one- 
by-one according to Wick's theorem. Even though this 
is not impossible from theoretical point of view (for ex- 
ample, the third cumulant has been derived using Wick's 
theorem^) , it is actually not practically executable. For- 
tunately, a successful solution to this problem has been 
given in Ref.H^ done by Gogolin and Komnik. In their 
seminal work)^^ they introduced a fictitious measuring 
field X in the tunneling Hamiltonian to count an electron 
when it goes through the system under studied in the 
direction of current and assumed it as a slowly- varying 
quantity in time (this is valid since it is set to be equal to 
zero in the final concrete calculations).^^ This presump- 
tion facilitates an adiabatic expansion of the generating 
function of the cumulants of charge current distribution 
to the first order. By doing so, the calculation of gen- 
erating function is then transferred to the calculation of 
a so-called adiabatic potential, which can be completely 
determined by the well defined single-particle Keldysh 
GFs. At the end, they obtained a generic formula for the 
generating function of the cumulants expressed only in 
terms of the local Keldysh GFs of the central region (and 
the tunnel-coupling and Fermi functions of the leads), as 
done in the Meir-Wingreen current formula?^ except for 
the presence of the counting field. The advantage of this 
Hamiltonian approach is that (1) it provides a system- 
atic and easy method to evaluate the cumulants of FCS 
by avoiding the tedious application of Wick's theorem; 
(2) more importantly, the derivation makes no assump- 
tions about interaction between electrons inside the cen- 
tral region; and (3) we can directly use the whole power 
of the Feynman diagram technique and connect to many 
known results of the NGF. Of course, a proper knowledge 
of the self energy is still indispensable for interacting sys- 
tems. With respect to this consideration, Gogolin and his 
coworkers first applied their approach to investigate the 
FCS of charge^SiiS current through a QD in the Kondo 
regime at the Toulouse limit, where the single- impurity 
Anderson Hamiltonian can be mapped to a quadratic 
form by performing a canonical transformation. More- 
over they analyzed the FCS of both charge^ and spin^ 
currents at the strong-coupling fixed point within the 
framework of the Nozieres-Fermi-liquid theory. Besides, 
they also studied the current cross-correlation (CC) cor- 
relations of a multiterminal Kondo-QD in the strong cou- 
pling limit. Another application of this approach is to 
explore the FCS of a molecular QD with strong electron- 
phonon interaction!^ 

However, it is known that the theories of the Toulous 
limit and the strong-coupling limit are both valid only 



at the deep Kondo region, where charge fluctuation is 
totally quenched and the dynamic properties of the sys- 
tem are completely determined by spin fluctuation. So 
in the present paper we will perform a complementary 
investigation of the FCS for a Kondo-QD in the multi- 
probe case at the deep and intermediate Kondo regions 
by means of the Coleman's infinite-?/ SBMFT ,-4^ which 
is believed to provide a proper description of the Kondo 
correlation at these regions even under nonequilibrium 
situation at zero temperature i^"— In combination of the 
SBMFT with the Gogolin and Komnik's approach, we 
find that the self-consistent equations, which determine 
the expectation values of the slave-boson fields, becomes 
counting field A-dependent. It is argued that these A- 
dependent terms of the mean values describe the Bose 
field fluctuation (it is equivalent to the charge fluctua- 
tion due to the completeness relation between fermions 
and bosons) induced by transport measurement. More 
interestingly, we find that these A-related parts of slave- 
boson fields generate an additional contribution in the 
zero-frequency shot noise formula, in contrast to our 
previous result without consideration of the Bose field 
fluctuations.'^^ Numerical calculations show that these 
additional term results in obvious enhancement of the 
current auto-correlation in the two-terminal setup and 
even a sign change of the current cross-correlation in the 
three-terminal setup of the Kondo dot in the intermedi- 
ate Kondo regions. 

The outline of the paper is as follows. In Sec. II, we 
give the model Hamiltonian describing electron tunnel- 
ing through an interacting QD attached to three leads in 
the presence of the three respective counting flelds and 
explain the theoretical formulation of the SBFMT, the 
counting-fleld-dependent self-consistent equations. In 
Sec. Ill, we provide the explicit expressions of the first- 
order differential equation of the adiabatic potential with 
respect to the counting fields, and the zero- frequency cur- 
rent auto- and cross-correlations at zero temperature. In 
particular, we discuss the additional terms due to fluc- 
tuation of slave-boson field. In Sec. IV, we perform con- 
crete numerical calculations and discussions for the cur- 
rent auto- and cross-correlations based on the formulae 
given in Sec. HI. Finally, our conclusions are given in 
section V. 



II. MODEL HAMILTONIAN AND 
THEORETICAL FORMULATION 

A. Model 

We model the electronic transport through a single- 
level QD coupled to multi-terminal using the infinite- [/ 
Anderson Hamiltonian as: 

■qka 
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rjka 

where c]^^^ {crfka) creates (destroys) an electron with mo- 
mentum k, spin a = {t,^}, and energy dispersion ejjka- in 
the lead ij, with 77 = {1, 2, 3} in the three-terminal case; 

(dcr) creates (destroys) a spin-cr electron on the QD 
with energy level e^; Jn is the coupling matrix element 
between the dot and lead rj; and A^(t) is the artificially 
introduced measuring field with respect to the lead rj on 
the Keldysh contour: A,,(t) — \-9{t)6{T ~t) on the for- 
ward path and A^(i) = \^9{t)0{T — t) on the backward 
path (7" is the measuring time during which the counting 
fields are non-zero and = — A,,+ = A)| iii40,42 

Under the framework of the infinite-?/ slave-boson 
approach, the ordinary electron operators on the QD 
are decomposed into a boson operator b (describing the 
empty state on the QD) and a fermion operator /o- (de- 
noting the singly occupied state with electron spin-cr), 
da ~ fab^ and dj. = In addition, a constraint, 

b^b + J2cr fa-fa- = Ij must bc imposed on these auxiliary 
operators as requirement of no double occupancy in the 
[/ — >■ cxD limit. Then, the effective Hamiltonian becomes 

r}ka (J 

-H^[7,e'^''«/V^6c,fc,+H.c.' 

rjk(7 

+^{b^ + Y.flf,-^, (2) 

with a Lagrange multiplier ^ to guarantee satisfaction of 
the constraint. 

B. Self-consistent Equations of Slave-Boson Mean 
Field Theory 

In the mean field approximation, the slave Bose opera- 
tors, w and 6, can be assumed as c-number and replaced 



where /») = [1 + exp (oj — ^.^i) j k-^T] ^ is the Fermi dis- 
tribution function at temperature T and chemical po- 



by their corresponding expectation values, h'^ and 6, in 
the Hamiltonian Eq. To determine the unknown pa- 
rameters, b and ^, we start from the constraint, and the 
equation of motion of the slave-boson operator using the 
Keldysh technique for systems out of equilibrium: 

Y,G-j + (t,t)+n, = 1, (3) 

(7 

e6 + E7.G7+,,,(i,i) = 0, (4) 

rfka 

where rn, — {wb) is the unoccupied-state number, 
GjJ'{t,t') = i{fl{t')fa{t)) is the dot non-equilibrium 

correlation GF, and 67+,,^^^ (i, t') = i{clka{t')fa{t)) de- 
notes the dot-lead correlation GF. With the Hamiltonian 
Eq. ([2|) , the mixed correlation GF can be readily cast in 
terms of Gjl^{t,t') with the help of the equation of mo- 
tion of the operators and then applying the Langreth an- 
alytical continuation rules in a complex time contour 4^ 
In this end, we obtain a self-consistent set of equations 
in terms of the QD's correlation GF in Fourier space: 

/ dioGj+{Lu)+n, = 1, (5) 
^ + ^ E / G7+(c.)(^ -6,-0-0. (6) 



Therefore, our next step is to calculate the dot cor- 
relation GF, GjJ'luj). It is evident that the mean-field 
Hamiltonian Eq. ([2]) is effectively a noninteracting reso- 
nant level model. By taking the counting field A,, {t) to be 
opposite constant on the forward and backward Keldysh 
branches as A,,_ — — A,,+ — Xr^^^ one can therefore read- 
ily evaluate the NGF of pseudo-fermion operator fa in 
terms of the original notation of Keldysh for GFs (the 
time-ordered GFs): 



(7) 



(8) 



tential (1^ — Ep + eV5, of lead i] [Ep is the Fermi en- 
ergy and Vr, is the bias- voltage applied to lead 77), and 



^f"^""' Vo{io){ -27n,V e-^^''/'r,(l-/,) -(c. - - ^ + V r,(2/, - 1) 



with (r = E,r,) 



Vo{oj) = (w - - + {ribTf + 4nl ^ r^r^,^(l - f^,) (e^(A.-A,0/2 _ 1^ 
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^r) = '"'Sfc ~ Eijfccr) is the coupling strength be- 

tween the QD and lead rj. In the wide band limit, we 
neglect the energy dependence of Tn and take it as a 
constant. Note that the GF formulae is similar to that of 
a noninteracting system, except with the effective energy 
level ed + S, and the effective tunnel-coupling constant 
ni,Tn instead, which renormalize the GF of the QD due 
to Kondo correlation under the approximation employed 
here. 

Substituting the resulting GF Eq. ^ into the self- 
consistent equations ([5]) and ([5]), we can obtain the two 
parameters, nj, and f , under a finite bias voltage. Obvi- 
ously, they are both functions of the counting fields A,,. 
For the sake of analysis and calculation of the shot noise, 
we can expand nt and ^ to the first order of as: 



(0) , • 
rib = nl ' + I 



n 

e = e(°)+^EJ^4'^+^(\1)' (10) 



where the superscripts (0) and (1) denotes the zeroth 
and first order terms of the expansion coefficients, re- 
spectively. The zeroth order terms, and ^^^^ are the 
original variational parameters of the SBMFT and irre- 
spective of the counting fields; while the first order terms. 



''bTj 



and are new here depicting fluctuations of the 
two parameters around their respective expectation val- 
ues due to measurement. 

The zeroth order self-consistent equations become 



.(0) 



duj 2(w - ed)(E^r').A) 



= 1, 



= 0, 



(11) 

(12) 



with id = £d + and F = T- 1* is noted that 
these two nonlinear equations are exactly the same as 
the previous results without the presence of the counting 
fields A^r^i^ and they constitute a closed set of equa- 
tions to completely determine the unknown parameters, 
n^°^ and for a particular QD system, e^, at a given 
bias voltage. The parameter characterizes the lo- 
cation of the Kondo peak in the quasiparticle density 
of states, and the parameter n[°^ mimics the width of 
the Kondo peak (i.e. the Kondo temperature Tk) and 
renormalizes the tunnel-coupling of the QD to the ex- 
ternal leads by Kondo correlation. It is known that the 
two parameters give the correct qualitative behavior of 
Kondo physics at zero temperature and low bias voltages 
{eV ^ Tji), and thus define the current-bias voltage char- 
acteristics and the differential conductance of a Kondo- 
QD] ^2i34,35,47 Furthermore, we can obtain the first order 



self-consistent equations as 

,(1) 



n 
)- 

n 

n 



(0) 
b 

(1) 



_ r(i)^(i) + r(i) - 

<^ (0) b CTj ~ '-'i 



where the integrals are as follows: 



1(0) ^ 



2vr (w-?,)2+r2' 



2TT [(^ - id)2 + r2]2 



t(i) - 



did 



2TT [(^ - g,)2 + r2]2 



2r») /jj 



(13) 
(14) 

(15a) 
(15b) 
(15c) 



27T [{u;- ?rf)2 + F2 

8('^r)'(E," r,./,.)E,, r,F,,(/, - /,,)] 



[{io - ?d)2 + F2 



(1) 



d^4[(..-6rf)2-F2](E^F,/,) 



27T [(^_g^)2+r2]2 



(15d) 



(15e) 



(15f) 



16(ni°))2(., 



duj 



ri" 



[(c.-?rf)2+F2]2 



(15^ 



These equations are one of the central results of this pa- 
per. Once n^^"* and are known, n^^^' and can be 
evaluated by solving the set of linear equations, Eqs. 
and (IT4)) . The two first order terms have no influence 
on the current but do affect the shot noise and of course 
the higher cumulants (see below). One can argue that 
the two terms depict the fluctuations of both the boson 
field b and the renormalization of the resonant peak ^, 
which can be ascribed to the charge fluctuation of the 
QD due to its attachment to external electrodes and ap- 
plication of bias voltage. Consequently, we will find in 
the following calculation that these terms have a nonzero 
contribution to the shot noise of the QD at a finite bias 
voltage. Moreover, they play a more important role at 
the intermediate Kondo regions, where the charge fluctu- 
ation is more profound, than at the deep Kondo regions. 
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III. FULL COUNTING STATISTICS BY 
SLAVE-BOSON MEAN FIELD THEORY 

In this section, we will first study the cumulant gener- 
ating function (CGF) x(A) = x(Ai, A2, A3) for the FCS of 
a Kondo-QD with three terminals within SBMFT. Then 
we will derive the formula for the current auto-correlation 
and cross-correlation taking into account the charge fluc- 
tuation effect in the two- and three-terminal configura- 
tions, respectively. 

The CGF can be calculated as a Keldysh partition 
functior 



,11 



x(A) = (Tce-^/c^^(*)^*), 



(16) 



where Tc is the Keldysh contour ordering operator and 



is the electron 



tunneling operator in the Hamiltonian Eq. ([2]). Accord- 
ing to Refs. Em to calculate the CGF x(A) it is tech- 



nically more convenient employing the adiabatic poten- 
tial method: Inx(A) — — iTZ^(A_,A+) = —iTU{\,—\), 
where the adiabatic potential Z^(A_,A+) is defined due 
to the nonequilibrium Feynman-Hellmann theorem as 



\ OXr,- I A 



4^ 

ka 



-/4 



7]ka 



H.c. 



(17) 



Therefore, the evaluation of the adiabatic potential 
amounts to a calculation of the mixed GFs, which can 
be cast into a combination of the dot GF G/o-(w) and 
bare lead GFs. By using the NGF obtained in last sec- 
tion and taking the counting field A,,_ = — A,;-|_ = A,j, the 
derivative of the adiabatic potential is given by 



dUiX,-X) 
dX„ 



27r 



I 



(18) 



where r,pj' (w) is the transmission coefficient of electron 
between leads 77 and ry': 



Trin' (w) = 



(19) 



It is clear to see that the transmission coefficient (1) has 
a simple Lorentzian line shape centered around the Fermi 
level (the Abrikosov-Suhl resonance) with a renormalized 
width due to Kondo correlation, and (2) is symmetric for 
exchanging 77 and 77': Tr,r,'{oj) = T,ffi{uj). It should be 
noted that due to the counting field dependence of the 
parameters Ub and ^ involved in the transmission coeffi- 



cient T^ri' (w), an explicit analytic expression for the CGF 
cannot be recovered here from Eq. (|18p . However, one 
can still obtain the m-th cumulant of the charge distri- 
bution involving lead r/ by executing the m-th derivative 
of the CGF, with respect to the corresponding counting 
field A,, at Ai = A2 = A3 = (we denote this condition 
as A = for shorthand in the following). In particular, 
the first two cumulants give the average current and the 
zero-frequency shot noise, respectively, which constitute 
our main objects in the present paper. For example, the 
current Ijj from lead ?/ to the QD is evaluated as follows: 



J 



^ ^ 2e 1 dlnxjX) 
^ h Td{iX^/2) 



4e dU{X, -A) 



dXr, 



A=0 



2e 
2e 



2tt 



A=0 



duj T,,,y {uj) {fn - fn' ) | a=o 
v' 



(20) 



It is clear that the three-terminal current formula Eq. (|20p is in perfect agreement with the previous result and 
reduces to the transparent form of the famous Meir-Wingreen current formula in the two-terminal case (for instance, 
by setting r3 = 0)41 
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We calculate now the zero- frequency current fluctuations of lead rj: 

4e2 1 d^\nx{\) 



5^,(0) 



h Td{iXr,l2f 



= 5(o)(0) + 5(^)(0), 



(21) 



A=0 



where 3),^} (0) is the mean-field result of the current correlation 

5^(0) = ^ fdJ J2 Tr,r,'{0j)[Ul - U) + /.'(I - /.)] 

and S^ij (0) is the correction to the mean-field result of the shot noise 



(22) 



A=0 



A=0 



(23) 



It is clear that sjpj{0) depends on the zcroth order terms of the two parameters, n^°^ and only, while "^^^^(O) is 



relevant to the first order terms, n[ll and and can be explicitly expressed as 



(w - ?d)2f^f^, 



n, 



(0) 



(a;-ed)2+f2 



2 



(24) 



The zero-frequency CC of currents through the different leads "q ^ -q' can be calculated as 

4e2 1 d'^\nx{\) 



Snn' (0) 



h T d{i\r,/2)d{iXr,'/2) 



= ^^°U0) + 5(;U0). 



(25) 



A=0 



Likewise, the CC £'^,,'(0) can be separated into two parts, the mean field result 5*^5(0) and its correction s'fl{Q): 



A=0 

(26) 



dT^m (^) 



d{i\,/2) 



ifr) fm ) 



A=0 



4e^ 



,(0) 



{uj - ed)2 + (f)2 



(1) (w--ed)IVV 



(a; - ear + (r)^ 



(/») /r/i)- 



(27) 



For illustration, we give the explicit expression for the mean- field shot noise, i.e., the current auto-correlation measured 
at the same lead (e.g., lead 1): 



^n'(O) = X y '^^{^i2[/i(l - /i) + /2(1 - /2)] + Ti3[/i(l - /i) + /3(1 - /a)] 

+Ti2(l - ri2)(/i - .f2f + Ti3(l - ri3)(/i - /a)' - 2ri2ri3(/i - /2)(/i - /a)} |a=o, (28) 
and the explicit expression for the mean-field current CC between leads 2 and 3: 

^23^(0) = "X / ^[^23[/2(l - /2) + /3(1 - fs)] + ^23(1 - T23)(/2 - hf + T^^T^ - h){h ' fl) 

+T2ir32(/2 - /l)(/3 - /2) + r3lT23(/3 - /l)(/2 - /3)] |a=0 • (29) 

I 



4e" 



Furthermore, in order to compare the present results with the previous formulae, we discuss the symmetrized 
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shot noise in the case of a QD connected to two reservoirs. 
By setting T3 = 0, the model employed above becomes a 
two-terminal system and the symmetrized shot noise is 
then defined as S{0) = [511(0) + 522(0) - 25'i2(0)]/4 
5(°)(0) + S'-'^^O), whose first part can be calculated via 
Eqs. ([22]) and ^ as 



S^'HO) = ^ J du{T,,[h{l ~ h) + f2{l - fi)] 

-TUh - f2f} ■ (30) 

Obviously, this formula is exactly the same as our previ- 
ous derivation applying Wick theorem within the frame- 
work of the This is the reason that we call it 
and Eqs. (12^ and (pS)) as mean-field terms of the current 
correlations. For reference, we briefly recall our previous 
derivation. The current operator flowing from the lead 
■q to the QD is dcflned as 



7]k<T 



, J2iVnb^cl,J, - V^nbc^ka), (31) 



r]k(7 



and the power spectrum of the CC correlator reads 



Sriri' ('^) 



dre— [({/,(t),V(i')})-(^)(V)], (32) 



with T = t — t' . Calculations of the correlator involves a 
pair of creation and destruction operators for each of the 
current operators at the two time variables. According 
to the conventional diagrammatic expansion technique of 
GF, the statistical expectation of the product of any two 
operators A and B can be divided into two parts 



{AB) - {A){B) + {AB),, 



(33) 



where (• • • )c means the connected part of the statistical 
average, which is vanishing for nonintcracting systems. 
As a result, under the framework of SBMFT and neglect- 
ing the fluctuations of both rtb and ^ (i.e. only consid- 
ering the leading terms, rS^^ and ^'•°''), we can eliminate 
the connected part of the above statistical expectation, 
and entirely contract the ensemble averaging with a prod- 
uct of four operators in ({/^(t), (i')}) as a product of 
ensemble averagings of two pairs of operators, such as 



iu.{t)u{t)cl,,At')Ut')) 



(34) 



and 



(35) 



All first terms in the left-hand side cancel out the term 
(/^ )(/,,') of Eq. (|32|) . whereas the second terms result in 
two lesser GFs running antiparallel in time (one starts at 
t going to t' , the other goes from t' to t), according to 
the Langreth theorem: the object with two antiparallel 
GFs C{t,t') — A{t' ,t)B{t,t') has the lesser component 
of C<(t,i') = A<{t',t)B>{t,t'). In succession, the hybird 
lesser GFs can further simplified as a proper combination 
of the local Keldysh GFs of the QD with the tunnel- 
coupling and Fermi functions of the leads, which leads 
in the end to the mean-field formula of the shot noise, 
Eq.®- 

Now we concentrate on the second part of the sym- 
metrized shot noise, Without loss of generality, 
we can assume Ai = — A2 (as a matter of fact, in the 
two-terminal case, although two counting fields Ai and 
A2 are introduced, the final result depends only on the 
difference A = Ai — A2 due to the charge conservation). 



and define 



,(1) - ^(1) 

'-62 



and ^{ 



(1) 



from Eqs. ([Ml) and we yield 



-e^'^ =^(^). Then 



4e^ 



8n 



(1) 



(w - e<j)2fif2 



{u - laY + (f)2 



(1). 



{uj - ed)TiT2 



ifi-f2). (36) 



From technical point of view, this term is a direct result 
of the connected part in Eq. (j33p . and consequently con- 
stitutes a correction to the mean-field result. Eg. (|30l) . 
Meanwhile, bearing in mind that in above mean-field 
derivation the fiuctuations of the boson field and the 
renormalization of the resonant level are both neglected, 
one can immediately conclude that this term 5'-'^^(0) is 
physically an account of the bias-voltage-induced fluc- 
tuation effect of the boson field. It is clear to see that 
S'^'^^O) is proportional to both n^^' and and is van- 
ishing under equilibrium situation. 



IV. RESULTS AND DISCUSSION 

In this section, we present the results of numerical cal- 
culations of the zero-frequency auto- and cross-current 
noises for the QD in the two- and three-terminal cases 
in the Kondo regime at zero temperature, paying special 
attention to the effect of boson held fluctuation on the 
current correlations, i.e., the correction terms of the shot 
noise formulae. The self-consistent parameters, n^j^^ and 



.(1) 



and can be obtained by solving Eqs. pT|) 



and ([T^ and ([Tl)) . respectively, for each bias volt- 

age and a given system e^. Using these calculated pa- 
rameters, we can then compute the mean-fleld results of 
shot noises and the total shot noises as functions of bias 
voltage according to Eqs. (HH), (UHl), and (gl]), re- 

spectively. In particular, we choose Fi = r2 = T/2 as a 



8 




(b) 



eV/T., 



Figure 1: (Color online) The calculated self-consistent pa- 
rameters, 71^^^ and nl^'' (a), ^'■"^ and ^'^^ (b) as functions of 
bias voltage V for the two-terminal QD with the energy level 
ea = -4.0. 



typical example and set F as the energy unit. The refer- 
ence energy will be always set at Ep = 0. For the mul- 
titerminal Kondo-QD, the width of the Kondo peak at 
the Fermi energy is characterized by the Kondo termper- 
ature, /cbTk = D exp (— 7r|e(i|/2F), at equilibrium {D is 
the energy cutoff and is set to D = lOOF here). In the 
following calculations, we confine the bias voltage from 
to SksTK/e {eV < SfcisTx) to ensure that the SBFMT is 
applicable for the Kondo correlation at nonequlibrium. 

First, we discuss the auto-current correlation for the 
two-terminal case, in which the external bias voltage is 
applied symmetrically, Vi = —V2 = V/2. For illustration, 
we depict in Fig. 1 the calculated self-consistent parame- 
ters, the zeroth order terms ^^'^^ and the first order 

terms nl^\ ^^^^ versus the bias voltage V for a QD with 
the discrete energy level = —4.0. It is easy to see that 
(1) these first order terms are nearly two orders of magni- 
tude smaller than their corresponding zeroth order terms 
at the whole region of the bias voltages; (2) they are in- 
creasing with increase of applied bias voltage; (3) they 
are vanishing at equilibrium, V — 0. These features are 
in good agreement with the consideration that the first 
order terms describe the fiuctuations of the correspond- 
ing parameters stemming from the bias- voltage-induced 
charge fiuctuations even in the Kondo regime. However, 
these small fiuctuations do cause an additional positive 
contribution to the mean-field shot noise as shown in 
Fig. 2(a,b), where the mean- field resuhs S'^^^O) and the 
total results 5(0) of the symmetrized shot noise are both 
plotted as functions of bias voltage for two QDs with 
Ed = —5.0, —4.0. Meanwhile, it is clear that the enhance- 



ment is more obvious in the system with = —4.0 than 
that of the system with ed — —5.0. This observation is 
consistent with the fact that the charge fiuctuation effect 
is entirely quenched in the deep Kondo regime but be- 
gins to play a unambiguous role in determining dynamic 
properties of the Kondo system when the system departs 
away from the deep Kondo regime. In Fig. 2(c,d), we plot 
the Fano factor 7 = 5(0)/2e/ [/ is the current calculated 
according to Eq. (|20|) ] as well. 




Figure 2: (Color online) The symmetrized zero-frequency shot 
noises (a,b) and the Fano factors (c,d) as functions of bias 
voltage in the two-terminal case for the Kondo-QD with sev- 
eral energy levels, ~ —5.0 (a,c) and —4.0 (b,d). The mean- 
field results and the total shot noises are plotted as solid-red 
lines and dash-blue lines, respectively. 

We now turn to the current CC in the case of three- 
terminal QD for Vi = -V2 = -V3 = V/2 and 
equal tunnel-couplings F3 ~ F/2. In Fig. 3, we show 
the zero-frequency CCs with and without the contribu- 
tion of the Bose field fiuctuation, 523(0) and ^23^(0), 
and the co rrespo nding Fano factors defined as 723 — 
823(0) /2e^/\hh\. It is well-known that in a noninteract- 
ing system, the sign of the current CC between different 
normal-metallic leads is always negative (antibunching) 
due to the Pauli exclusion principle of Fermionic statistics 
of electrons?^ which has been confirmed experimentally 
in a Hanbury-Brown-Twiss setup (HBT).— Nevertheless, 
a positive current CC (bunching) has been predicted 
in certain situations: a hybrid superconductor-normal 
system!^ spin-dependent sequential tunneling through a 
QD with three ferromagnetic leadsf^ a Coulomb interac- 
tion coupled systemi^i^ Very recently, a sign change of 
the current CC due to dynamical channel blockade has 
been experimentally observed in a subtle measurement of 
the CC between two output terminals in sequential tun- 
neling through two capacitively coupled QDs connected 
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to four independent electrodes 



54 




Figure 3: (Color online) The current CCs between leads 2 and 
3, 5*23(0) (a,b), and the corresponding Fano factors 723 (c,d) 
as functions of bias voltage in the injecting lead in the two- 
terminal case for the Kondo-QD with several energy levels, 
Ed = -6.0 (a,d), -5.0 (b,e), and -4.0 (c,f). The mean-field 
results and the total shot noises are plotted as solid-red lines 
and dash-blue lines, respectively. 

For the case of Kondo-correlated transport, we notice 
that there have been two theoretical work devoted to 
the current CC in this regime and no sign reversal has 
been found.— Schmidt and his cowork have investi- 
gated the FCS of Kondo-type tunneling through a multi- 
terminal QD using the Nozieres-Fermi-liquid theory and 
reported analytical expressions for the current CCs at the 
true strong-coupling fixed point and the departure from 
this point 45- Sanchez and Lopez has applied the same 
SBMFT as the present work to calculate the current CC 
of a three-terminal QD in the presence of ferromagnetic 
leads.-** Different from this paper is that they take no 
account of the fluctuations of both the Bose field and the 
level renormalization. In other words, their calculations 
correspond to the mean-field results of the present inves- 
tigation. Therefore, our mean-field Fano factors 723 as 
shown in Fig. 3(d-f) by the solid-red lines has the same 
properties as those in Ref. [sj in the case of paramagnetic 
leads: (1) it is always negative; (2) it has a minimum at 
V = and the value of the minimum is nearly equal 
to —0.44; (3) it reaches a saturation value at very high 
bias voltage. It is however worth noticing that a rate 
equation calculation has addressed a sign change of the 
current CC due to dynamical spin blockade for the se- 
quential tunneling through the same setup, i.e. a QD 



connected to three ferromagnetic leads<^ What is inter- 
esting in the present calculations is that we indeed find 
a positive current CC for the QD with = —4.0 under 
higher bias voltages V > 2.0Tk, as shown in Fig. 3(c,f) by 
the dash-blue line, if the charge fiuctuation is taken into 
account, even though the system under studied involves 
no ferromagnetic electrode. But, the current CC remains 
negative for the systems with = —6.0 and —5.0, since 
the charge fluctuation is weaker in the two systems than 
that in the system with = —4.0. Generally, the overall 
role of charge fluctuation is to reduce the fermionic HBT 
type correlation at higher bias voltage region, which is 
consistent with the analytical conclusion in Ref. 



V. CONCLUSIONS 

In conclusion, we have analyzed the effects of charge 
fluctuation on the current statistics in a Kondo QD con- 
nected to multi-leads at zero temperature on the basis of 
the infinite-?/ SBMFT to demonstrate the Kondo cor- 
relation. By introducing counting fields with respect 
to each electrodes and employing the NGF technique, 
we have derived an explicit analytical expression for the 
first-order differential equation of the adiabatic potential 
with respect to the counting fields. Meanwhile, we have 
found that the SBMFT self-consistent equations become 
dependent on the counting fields, which provides a pow- 
erful instrument for taking account of the effect of the 
Bose field fluctuation, i.e. the charge fluctuation effect, 
at nonequilibrium on the current correlation. In particu- 
lar, we have performed derivations for the concrete ana- 
lytical expressions for the symmetrized shot noise in the 
case of two-terminal QD and the current CC in the case 
of three-terminal QD. It has been found that the con- 
sideration of the effect of charge fluctuation generates an 
additional contribution to the mean-field result, which 
leads to enhancement of the symmetrized shot noise in 
comparison to the mean-field result and a positive cur- 
rent CC for certain system parameters. 
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